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ON THE DEFINITION OF CRATER SATURATION DIAMETER 


BY 


Allan H. Marcus 


ABSTRACT 


A theoretical model is developed for determining the 


diameter D such that the density of lunar craters exceeding 
P 

constitutes a specified proportion p of the maximum number 
that could be observed on a crater- saturated surface. It is 


shown that this diameter is sensitive to the value of the exponent 


in the fresh crater size distribution, and is much more informative 
than the usual log-log plot intercept definition of saturation 


diameter. 
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ON THE DEFINITION OF CRATER SATURATION DIAMETER 

by 
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Theory 

In order to clarify the definition of a planetary surface 
which is "saturated" by craters, we must consider unsaturated surfaces 
as well. We start therefore with the expected number of craters of 
diameter D to D + dD formed per unit area per unit time on a 
surface, n o (D)dD. If the crater centers are distributed at random 
on the surface, and if the formation of a crater of diameter y will 
destroy a crater of diameter D if the center of the y-crater falls 
within an area A(D,y) centered on the center of the D-crater, then the 
net average destruction rate of D-craters is A(D), 


A(D) = f A(D,y) n (y)dy 0-) 

J 0 ° 

Hence the expected number n(D,t)dD of craters of diameter D to 
D + dD, per unit area per unit time at time t, satisfies the 
differential equation 

~ n(D,t) = n (D) - A(D)n(D, t ) (2) 

at o 

and with the initial condition 


n(D,0) « 0 



2 


we obtain 


n(D,t) = n(D, a>) [1 - exp(- A(D)t)3 


where the equilibrium number density is 


n(D,=) = n Q (D) / A{d) 


It can be shown (Marcus, 1970) that under a very broad range 
of crater obliteration models, there is an equilibrium density C 
such that if for some s > 2, D > 0, cp > 0 , 


then 


where 


n (D) = scoD 3 D" 3 " 1 
o ' o 


A(D) = k if 


n(D, <») = 2 C D -3 


C = scp D q /2k. 


Eq. (3) can be recast into the form of a cumulative expected number 
density N(D,t), using (4)-(8): 


N(D,t) = f n(D,t)dD 


t(b+i 


f 1 ** -< «j ] 
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where 

40 a 

r{B+l) = J u p exp(-u)du 
0 

6 = 2/ (s-2 ) 

e = kt v>~ 2 ^ , D = (kt. / € ) 

p{ Y < (r] = f u^ _I exp(- u )du/r(3) 

p “ '0 

Hence: A SURFACE is lOOp $ SATURATED AT DIAMETER D p IF 

N(D pS T) = p M(D p ,«) 

i.e. 

(i - p ) ^ / r(P + i) - pf \ p < c p } (i°) 

Gault (1970) has used another definition. Some numerical examples 
are instructive. 

EXAMPLE 1. p ~ 0 . 

e must be small, so expanding the exponential in 

' P 

pfv < <7 } as a MacLaurin series in u, we obtain 

1 T 6 - - p J 

e p „ (p + l) p/6 < u > 


D p „ (pkt/(p + l)p ) 9//2 = C p 


( 12 ) 
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EXAMPLE 2. The usual definition of the diameter D at which a 
surface is saturated (Gault et. al. , 1970) i- s 


N(D , <» ) = t f n (D)dD 


J * 0 


whence, from (5) and (7 ) . 


* 

D = C n 


where C is defined by (12) with p = 1. However, it is evident 


from (3) that no surface is ever 100$ saturated. 


EXAMPLE 3- p ~ 1 

must be large, so <sj ~ 1 and 


e p ~ (r(B + l) / (l - p)) 


~ (kt) 0 / 2 ((i - p)/r(B + 1)) 1//2 = 


EXAMPLE k. As a preliminary estimate of D , use saturation ratio 

i: 

for the incremental number density 


p = n(A p ,t) / n(A p ,») 


whence from (3) a nd (6) 


k - ( - kt/log (l - p)) 
P ^ 


(18) 
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We obtain after some numerical calculations 


€ p 

S - k 

s = 3 

i s = 8/3 

p . 1 

8=1 

8 = 2 

1 8 = 3 

0.5 j 

2.25 

| 1.15 

1 

i 0.97 
\ 

0.8 j 

4.965 

! 2.73 

1 2.25 
» 


Thus : 


p 1 

s = 4 s 
8 = 1! 

s = 3 ! 
8 = 2 

s = 8/3 = 2.67 
8 = 3 

C /D • 
p p ' 

0 . 5 j 

H 

\n 

0 

1.53 ; 

1.75 | 

i 

0.8} 

1.74 j 

t 

2.28 1 

3.06 : 

- — 1 

i 

* / ' 
D /D i 
P ■ 

0.5 

1 

1.06 | 

0.767 ! 

i 

.... . 

i 

0.620 | 

1 

j 

! 

j 

0.8 

1.55 

1.82 

2.18 ! 

} j. 

B /D 1 
D P P 

0.5 i 

1.06 

1.15 

i i - 6 5 ; 

0.8 

i 

j 0.996 

j 0.862 

| 0.616 

A /D 
P P 

0.5} 

1.80 

1.66 

| 1.65 1 

> * 


0.8 

; 1.75 

1.70 

i 1.65 

i 

? 

1 


Discussion 

The results are rather sensitive to the choice of a 

definition for crater saturation. D , the diameter such that the 

number of large craters exceeding this size is 100 p percent of the 

total number that could be observed, is a well-defined, computable, 

and intuitively satisfactory characterization of saturation. The same 

could be said of the diameter A defined by incremental number 

P 

densities, and it can be shown that if s not too much larger 
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than 2, then 

A p _ (log e 2) D p - 1.61 D p 


(19) 


as the numerical example shows. 

Both of these quantities must be estimated from an empirical 
analysis of crater densities, based on the result 


N(D,t) ~ C/D 2 


if D < < D 

P 


2Ckt / sD S 


if D > > D 

P 


( 20 ) 


We usually plot log W(D,t) vs. log D and define D as the 
intersection of these two (approximately straight line segments. 
The curved transition between the straight line segments covers a 
large range of diameters and is obscured by the logarithmic plot. 


It is interesting to compare crater counts from the 
Apollo 11 and 12 crater count (Shoemaker, 1970), for D in meters; 


i 

s i 

C 1 

2Ckt/s : 

D* 1 
i 

3 

D 0.5 i 

D 0.8 

Apollo 11 j 

2.93 | 

0.!0/ m 2 j 

8.0/ 2 
m 

: ikl 

m i 

2.15 

l8 9 m 

75 m 

( 

Apollo 12 

! 2,86 ! 

i • 

0 . 10/ b 2 ] 

2-0 / 2 
1 m 

32.6 1 
• m 

i 

2.33 

4 

45 - 5 m 

| 16.8 

1 m 

j 


The D criterion is evidently more informative. 
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